Experimentally Realizable C-NOT Gate in a Flux Qubit/Resonator System 
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In this paper we present an experimentally realizable microwave pulse sequence that effects a 
Controlled NOT (C-NOT) gate operation on a Josephson junction-based flux-qubit/resonator sys- 
tem with high fidelity in the end state. We obtained a C-NOT gate process fidelity of 0.988 (0.980) 
for a two (three) qubit/resonator system under ideal conditions, and a fidelity of 0.903 for a two 
qubit/resonator system under the best, currently achieved, experimental conditions. In both cases, 
we found that "qubit leakage" to higher levels of the resonator causes a majority of the loss of 
fidelity, and that such leakage becomes more pronounced as decoherence effects increase. 



PACS numbers: 74.50.-|-r, 03.67.Lx, 85.25.Cp 



I. INTRODUCTION 

The basic requirements for a successful quantum com- 
puter have been expressed by David DiVincenzo nearly a 
decade ago.^^ These five criteria, listed below, have been 
widely accepted as being the best road map for achiev- 
ing realizable quantum computing by most research pro- 
grams throughout the world: 

1. A scalable physical system of well-characterized 
qubit s. 

2. The ability to initialize the state of the qubits to a 
simple fiducial state. 

3. Long (relative) decoherence times, much longer 
than the gate-operation time. 

4. A universal set of quantum gates. 

5. A qubit-spccific measurement capability. 

At this time, there are various schemes being proposed 
to satisfy the above criteria and realize a quantum com- 
puter4 At the few-qubit level, these schemes include 
those based on trapped ions;^ liner opticS)^!^ and nu- 
clear spins in liquid-state moleculesiii^ For the long-term 
prospects of scalability though, those schemes that utilize 
Josephson junction-based qubits^ have significant advan- 
tages with current solid-state manufacturing technology. 

Since the initial breakthrough in the coherent manipu- 
lation of a single Josephson junction-based charge qubit 
nearly a decade agO;^ the experimental focus has mostly 
shifted to the creation, control, and subsequent manipu- 
lation of, multi-qubit entanglement in similar Josephson 
junction-based systems. For example, coherent oscilla- 
tions between two qubits have been observed by using a 
fixed inter-qubit coupling fiSiiiii^ However, the fixed na- 
ture of the qubit-qubit coupling used in these experi- 
ments makes it difficult to scale up such circuits in the 
future. To overcome this problem, a fast switchable cou- 
pling between two qubits has been propose d^^i^^ and also 
demonstrated^^ 



Beyond direct qubit-qubit couplings, another solution 
is to make use of a quantum bus (qubus) as a coupler 
between qubitsj ^^i^^'^^ Using the "qubus" concept, we 
can perform any two qubit operation, between any two 
qubits that are coupled to the bus, without using multi- 
ple swap gates, which are necessary in other systems us- 
ing direct qubit-qubit coupling. In particular, harmonic 
oscillators formed by superconducting circuits seem to 
be a good candidate for a "qubus" -type coupler. Early 
experiments with such couplers have shown coherent os- 
cillations between a qubit and the oscillator, which was 
made of lumped elements, namely capacitors and in- 
ductorsJ^i^ More recently, a distributed circuit, based 
on a coplanar waveguide resonator, attracted consider- 
able attention as an oscillator because of its high quality 
factor (Q-factor) and impedance matching to other cir- 
cuits i^iS^ Furthermore, coherent quantum state transfer, 
between two Josephson junction-based qubits, via such 
a waveguide resonator, has been demonstrated in both 
the phase^^ and charge regimci^i Lastly, more recently, 
two-qubit algorithms have been demonstrated in a two- 
transmon qubit/resonator systemi^ 

Because of the experimental viability of the "qubus" 
coupler concept, as well as its obvious advantages in scal- 
ability, we are using this paper to propose an experimen- 
tally realizable microwave pulse sequence that will enact 
a controlled NOT (C-NOT) gate between two supercon- 
ducting flux qubits, which are coupled via a harmonic 
oscillator bus. This pulse sequence flux qubit/resonator 
design is complementary to other recently proposed sys- 
temSfi^iS^ but stands apart due to its efficiency of oper- 
ation and fidelity of its end state. 



II. EXPERIMENTAL SYSTEM 

One of the most promising solid state quantum com- 
puting elements is a superconducting flux qubit that typ- 
ically consists of three Josephson junctions in a loop: 
Two of equal size, one smaller by a factor a ~ O.Si^l 
The sizes of the junctions are chosen so that the goo- 
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metric self-inductance of the loop is not physically rele- 
vant. The two lowest energy states of the qubit at the 
flux degeneracy point, usually denoted as |0) and |1), arc 
superpositions of macroscopically distinct clockwise and 
counter-clockwise persistent current states. Repeated ex- 
periments over the past few years have shown that this 
type of flux qubit is a well-deflned quantum system that 
can perform singlc-qubit rotations22i22ii2£i2ii22 as well as 
achieve longer coherence times than other kinds of su- 
perconducting qubits.— However, we have to operate the 
flux qubit at the degeneracy point in order to enjoy the 
best coherence time because it is easily affected by a flux 
fluctuation apart from the degeneracy point i^i^i^^ 
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FIG. 1: (Color online) Multiple flux qubits (blue rectangles 
with three crosses), each addressed by its own microwave line 
(green line), coupled to a resonator via a mutual inductance 
M. The resonator is schematically represented by the induc- 
tor L and the capacitor C. 

With the above in mind, we are currently working with 
the following architecture (see figure [1]), built around the 
previously discussed three Josephson junction qubit de- 
sign, which incorporates the fundamentals of the "qubus" 
concept Each qubit couples to a resonator through a 
mutual inductance M. Here the resonator is schemat- 
ically represented by lumped elements, but it can be a 
distributed circuit, for example, a coplanar strip waveg- 
uide. External magnetic flux through each qubit is a half 
flux quanta to set the qubits at the degeneracy point. 



III. C-NOT PULSE SEQUENCE 

In order to execute a C-NOT gate in our architecture 
(see figure [T]) we have to entangle the flux qubits with the 
resonator. To do that, we need to apply a sequence of 
DC-shift pulses, or apply a sequence of microwave pulses, 
through microwave lines to our qubits. In the case of 
the former pulses, we can adiabatically change the qubit 
frequency to fit the resonator frequency, but also non- 
adiabatically to create a coupling between the qubit and 
the resonator. As a result, we can "turn on" the coupling 
between the qubit and the resonator non-adiabatically, 
making them into an entangled state j^°i^^ However, these 
pulses cause DC-based excursions away from the flux de- 



generacy point, and can reduce the dephasing times of 
the qubits drastically. The large bandwidth of the pulses 
can also reduce the overall gate fldelity. Fortunately, a 
special flux qubit design may solve these problem by us- 
ing elaborate DC pulsesi^ 

In the case of the latter pulses, we can create entan- 
glement between the qubits and the resonator by us- 
ing a known two-photon blue sideband (BSB) transition 
at the qubit 's degeneracy point J^i^^ These microwave 
pulses have a more narrower bandwidth than the DC- 
shift pulses, thus allowing us to obtaining a much higher 
fidelity gate, as well as minimizing pulse-induced dephas- 
ing. C-NOT gate operations based on such BSB tran- 
sitions have been achieved with a high fldelity in ion 
trap experimentsi^i^ In these experiments, an elaborate 
controlled- gate, using four BSB pulses, was used4 

It is an interesting idea to introduce such a BSB- 
based C-NOT gate operation to our superconducting 
qubit/rcsonator system. However, it is hard to achieve 
a similar high fidelity because of the strong fixed cou- 
pling between the qubit and the resonator. This coupling 
leads to larger energy shifts at higher energy levels [see 
figure[2ta)]. Indeed, the BSB-based C-NOT gate uses up 
to n = 2 states, where n is the photon number in the res- 
onator. Therefore, we will utilize this coupling to realize 
a high-fidelity controUed-Z gate instead by only using up 
to n = 1 states. From this, we can build a pulse sequence 
that wiU do a C-NOT gate [figure [DJb)]. 

As figure [DJb) shows, the first BSB tt pulse transfers 
information from the control qubit to the resonator, and 
the second one transfers the information back to the 
control qubit. During a free evolution between the two 
pulses, we use only four energy levels, |10) and |11), for 
n = and n = 1, where \CT) = \C)\T) and |C) {\T)) 
represents the state of the control (target) qubit. The 
energy difference of the target qubit between n = and 
n = 1 indicated by the two dashed arrows in figure [2Ja) 
enable us to realize a controUed-Z gate between the res- 
onator and the target qubit with a suitable temporal de- 
lay between the two pulses. Similar gate operations have 
been employed in nuclear magnetic resonance quantum 
computationi^ These pulses, with proper phases, work 
as a controUed-Z gate between the two qubits (C/cz)- 
We have to adjust the phase difference between the two 
pulses carefully to cancel the effect of the off-resonant 
ac-Stark shift S in figurel^Ja) Finally, four single-qubit 
gates, and the Ucz gate, form the C-NOT gate (C/c-not) 
as follows: 



NOT 



10 
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(1) 



in the two-qubit basis {[11), |10), |01), |00)}. Here, Yq 
and It represent a ^ rotation of the control qubit, and 
that of the target one, respectively. 
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(a) 




FIG. 2: (Color online) (a) Energy levels of the two 
qubit/resonator system and transitions used in our pulse se- 
quence. Solid (dashed) lines represent the energy levels with 
(without) the coupling between the qubits and the resonator. 
Additional shifts 5 occur during the irradiation of the BSB 
pulse. Energy difference between the dashed arrows is utilized 
to perform the controUed-Z gate between the resonator and 
the target qubit. (b) Basic operational sequence of "target," 
"control," and BSB pulses on our two qubits, using values 
defined in section HV Al that enact a C-NOT gate in our sys- 
tem. Each pulse is represented by the microwave frequency 
(Carrier frequency or BSB frequency), the rotation angle, and 
the phase of the microwave pulse. Note that the parameters 
T, the free evolution time, and <j), the microwave pulse phase, 
are variable: They are defined by the experimental system to 
those values that optimize the operation of the gate. 



IV. NUMERICAL SIMULATIONS 

As has been noted recently;^ in order to accomplish 
our C-NOT, or any usable quantum computational gate, 
it will be necessary to manipulate multiple numbers of 
two-level quantum states, which will be our physical 
qubits, at will with some type of pulse sequence. For our 
system, as discussed in scction lllll and schematically rep- 
resented in figure [H this means that any single- or multi- 
qubit gate operation will be done using tt pulses, ^ pulses, 
and so on, most likely delivered by an on-chip microwave 
line, or lines. Our architecture not only achieves this, but 
also allows for a Hamiltonian representation that easily 
allows us to simulate whatever necessary pulse sequences 
we see fit to enact. In this section, our focus will be on 
simulating those pulses that implement our high-fidelity 
C-NOT gate upon any two qubits in the architecture. 

To begin, when all the qubits in our system are at 
their respective degeneracy points, the overall system is 



represented by the following Hamiltonian: 

H = hfycsi^a + '^{^hfkCTzM + hgkax,k{a) + a) (2) 

k 

+hAMW,k sin(a;fet 4- (t)k)(Tx,k)- 

Here, h is the Plank constant, /res is the frequency 
of our resonator, \hfk(Tz,k represents our "k-th" qubit, 
hgka'x,k{0'^ + a,) is the coupling term between the "k-th" 
qubit and the resonator, and hAMW,k sin([jfci -|- (j)k)<^x,k 
describes the microwave pulse for each qubit. Since we 
are only looking at C-NOT gates, we will define one qubit 
as the "target," and one as the "control" qubit. Lastly, 
we should note that, in principle, our gate operation uses 
only the ground, and first excited state, of the resonator. 

Now, because we are wanting to analyze the time- 
evolution of multiple qubits - in particular, whether or 
not our pulse sequence accurately describes a C-NOT 
process on any two of them - we naturally will want 
to work within the density matrix formalism. In more 
detail, our simulated two-qubit gate, including error 
sources, is properly described by a completely positive 
map e, the density matrix of the output state pout of 
which can be written in the operator sum representation 
as^ 

Pout = e(Pin) = ^ E„iPinElxmn, (3) 

where Em are operators forming a basis in the space of 
4x4 matrices and is the density matrix of the input 
state. This expression shows that e can be completely 
described by a complex number matrix, Xj once the set 
of operators has been fixed. 

In this work, we have chosen these 16 operators to be 

Eii+j =A^(E) Aj (4) 

in the case of a two-qubit gate. Here, ~ I, Ai = 
f^x, A2 = (Ty, and A3 = Oz- To determine x, which is 
a 16 X 16 matrix, we need to simulate our gate for 16 
linearly independent input states IV'in)- For these states 
we have chosen 

IV-in) = ® IV-j), (5) 

|V,)e{|0),|l),(|0) + |l))/V2,(|0)+z|l))/^/2} 

as the initial states. From this we can "measure" the 
gate fidelity using the process fidelity 

■Pp = Tr(xidcalXsim), (6) 

where Xideai and Xsim represent the ideal matrix, and the 
one obtained from the simulation, respectively. 

A. Simulation Parameters 

Our two qubits (the "target" and "control") and res- 
onator operating frequencies, as well as couplings, were 
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defined in the following way: 

/l = /control = 6 GHz 
/2 = /Target = 5 GHz 

/,es = 10 GHz. 



gi = 0.1 GHz 
32 = 0.1 GHz 



(7) 



Wc tlien generated numerical simulations of equation ([3]) 
using the following values for our target and control 
pulses: 



Target Pulse 
Amw.2 = 0.1 GHz 
W2/27r = 4.99875 GHz 
(i = or TT 



Control Pulse 
^MW4 =0.1 GHz 
cji/27r = 5.9981 GHz 
= 0, 



(8) 



(9) 



as well as our blue-side-band (BSB) pulses: 

Control BSB 
^Mw,i = 2 GHz 
uji/2t: = 7.5601 GHz 
fli = or 0.347r. 



Here the rise/fall time of each pulse was set at 0.8 ns 
and the duration of the carrier 7r(^)-pulse was set at 
5(2.5) ns. As we mentioned in section Hill we optimized 
the duration of the BSB pulse, the free evolution time 
T, and the phase of the second BSB pulse in order 
to achieve the best gate fidelity; obtaining 16.295 ns, 
162.865 ns and 0.347r, respectively. The total length of 
the pulse sequence is about 200 ns, which is much shorter 
than the reported coherence time of a flux qubiti^ Lastly, 
we assumed that the first five resonator levels were acces- 
sible, thus n = to n = 4, where n is the photon number 
in the resonator. 



V. OBSERVATION AND ANALYSIS 

Multiple simulations, using the previously defined pa- 
rameters with various pulse shapes and operation times, 
have yielded an experimentally realizable matrix of tar- 
get, control and BSB pulse sequences, represented in 
figure El^b) , that will initiate a C-NOT gate upon our 
two qubit/resonator system with high process fidelity, 
Fp — 0.988, and within the coherence time of known 
flux qubit systems. The computed gate is as follows (to 
four significant figures): 



^c-not 



/ 0.0001 0.9863 0.0013 0.0022' 
0.9863 0.0001 0.0020 0.0013 
0.0021 0.0014 0.9886 0.0034 

V 0.0013 0.0023 0.0035 0.9880. 



(10) 



We also evaluated the effects of decoherence on our C- 
NOT gate by introducing a linear loss to the resonator 
(quality factor Q), as well as relaxation rates Fi and de- 
phasing rates r2 to the qubits, via a master equation of 



the Lindblad form (see figures [3] and [4]). Here we as- 
sumed that Pi and r2 for both qubits were equal. We 
obtained a process fidelity of 0.903 in the best conditions 
which have been achieved experimentally, for example, 
Q = 106, Pi = Pj = 0.25 MHz. Although the fidelity 
without the decoherence is not unity, it is much better 
than that with the decoherence, hence our pulse sequence 
could be useful enough to demonstrate the C-NOT gate 
exp er iment ally. 
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FIG. 3: (Color online) Process fidelity as a function of qubit 
decay rate. The upper (lower) curve represents the fidelity 
when the quality factor of the resonator is 10^ (3 x 10^). 

Now, in the case of flux qubits, it is difficult to fabri- 
cate a qubit which has an exact designed gap frequency at 
the degeneracy point, unless using the qubit with a con- 
trollable third junction.— Hence, we also simulated the 
case in which the two qubit frequencies became closer to 
each other. When /ics = 10 GHz and /control = 6 GHz 
and /Target = 5.5 GHz, we obtained a process fidelity of 
0.986 without decoherence. Decreasing the difference be- 
tween two qubit frequencies did not significantly affect 
the fidelity. 

This gate also works with any two qubits in a sys- 
tem with many qubits, since we can decouple any un- 
used qubits from the resonator via very quick (« 5 ns) tt 
pulses applied at the center of the two BSB pulses. For 
example, we obtained an process fidelity of Fp = 0.980 
in a three qubit/resonator system where the frequency 
of the third unused qubit was 7 GHz, its coupling to the 
resonator was 0.1 GHz, and the microwave frequency of 
the decoupling pulse was 6.9973 GHz. The other sys- 
tem conditions were similar to that use in the two-qubit 
simulation. 

Unfortunately, our gate fidelity, while still better than 
that currently achieved by other experimental systems 
(see for example Riebe et. al. and their table II) is still 
not unity. The lack of unit fidelity in our C-NOT gate is 
due to "qubit leakage"; the unintended, and persistent, 
excitation of higher-order resonator levels. As Fazio et. 
al. noted in 2000 qubit leakage issues are due to the 
fact that the qubit Hilbert space is actually a subspace 
of a larger Hilbert space; in our ease, one that includes 
accessible higher-order resonator levels. Although their 
work focused on charge qubits, their basic conclusions 
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hold for flux qubits as well, as our simulations have shown 
(see figure HI). 





protection through non-demolition detection and state 
re-injection^^i^i^i^ correction by teleporting "leaked" 
qubits back into the 2-level subspace^^ and specialized 
leakage reduction units (LRU's) which suppress leakage 
through a pulse sequences in a similar manner to bang- 
bang protection for qubit dephasingj ^^i^^ Unfortunately, 
these techniques arc quite cumbersome and involve spe- 
cialized protocols that go beyond the standard operations 
to realize QEC. 

Ideally, for large scale qubit applications, we do not 
want to correct for qubit leakage, instead we wish to sup- 
press it to a sufficient level where it can be ignored. It 
should be noted that this requirement is not as trivial 
as it might appear. For large scale applications, fault- 
tolerant quantum computation requires extremely low 
error rates. This "target" error rate for fault-tolerant 
QIP, often referred to as the "fault-tolerant threshold" is 
highly dependent on the underlying physical architecture 
and can be anywhere between 10"^ and io~7 ^47,48 rpj^g 
important point to realize is that if we wish to ignore co- 
herent leakage errors, and avoid implementing specialized 
protocols, errors caused by qubit leakage will need to be 
several orders of magnitude below any other correctable 
source of qubit error. Therefore, the pulse design of this 
bus mediated qubit/qubit coupling will be required to 
exhibit effectively zero leakage, even before open system 
quantum effects are even considered, if it is to be used to 
do realizable QIP. This means that a new pulse sequence 
will need to be created. 



FIG. 4: From top: Initial (left) and final (middle and right) 
density matrices (real components only) showing the |00) i-^ 
|00), |01) ^ |01), [10> ^ 111), and |11) ^ |10> operations, 
and the corresponding qubit leakage into higher resonator 
modes (n > 0), due to the pulse sequence values defined 
in section FlV Al without and with decoherence [Q — 10* , 
Fi = r2 = 0.25 MHz). The color scheme used is a nonlinear 
gradient: Zero values are white, with negative values greenish 
and positive values reddish. 

In more detail, coherent leakage into higher order 
modes of the harmonic oscillator bus seems to be respon- 
sible for the missing process fidelity of the system. Not 
only is this large an error rate unacceptable from the 
standpoint of large scale, fault-tolerant quantum com- 
putation, but the error channel itself is problematic to 
correct under current quantum error correction (QEC) 
schemes. 

One of the fundamental assumptions of QEC is that 
errors affecting qubits are localized to the two-levels as- 
sociated with the qubit encoding. Once errors begin to 
violate this assumption, either through state leakage or 
actual qubit loss, additional mechanisms must be em- 
ployed, at the encoding level, to allow for correction. 
Many different techniques have been developed to com- 
pensate for this type of error channel, including leakage 



VI. FOLLOW-ON PULSE SEQUENCE 
POSSIBILITIES 

From our work, it is our belief that the seen qubit 
leakage, even taking our simple decoherence model into 
account, and the corresponding non-unity fidelity is an 
unavoidable consequence of the programmed waveforms, 
and expected physical implementation, of our microwave 
pulses. Experimentally, in order to avoid qubit state 
transitions, the rise and fall times of our pulses have 
to be long enough to fulfill the adiabatic condition, but 
also short enough to avoid unwanted relaxation processes. 
Temporal lengths of the various target, control and BSB 
pulses used also depend on the microwave pulse ampli- 
tude, and are experimentally determined from the sys- 
tem's driven Rabi oscillations. Thus, while these numer- 
ical simulations demonstrate the principal application of 
a bus- mediated controlled qubit/qubit interaction, the 
implementation of the proposed pulse sequence requires 
further improvement of the leakage to build up not only 
a several qubit system, from the standpoint of a quantum 
testbed, but also a really-scalable QIP system. 

In more detail, analysis of our control, target and BSB 
pulse matrices indicates that our qubit leakage is initi- 
ated by our first BSB pulse on the control qubit. Since 
the second BSB pulse on the same qubit is only different 
by a phase, we only require a modified version of the BSB 
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pulse matrix to avoid leaking population into higher or- 
der modes of the harmonic oscillator. Also, as the simple 
"top hat" pulse design we have examined is already very 
well confined, it is expected that this pulse represents a 
good first approximation to a zero leakage modification, 
and due to the structure of this sequence, we already have 
a good idea which section of the pulse is "leaky." Thus 
we feci that the only way to truly generate a microwave 
pulse-induced gate, with fidelity experimentally equiva- 
lent to unity, for this system will be to use technologies 
based on optimal quantum control 

Optimal quantum control is a vast area of theoretical 
and numerical techniques that were developed largely in 
the field of nuclear magnetic resonance (NMR) and quan- 
tum chemistry, in order to find complicated control fields 
for NMR manipulations, or to increase chemical yields 
in chemical reactions involving photonic reagents. This 
research has more recently been applied to the quan- 
tum information field j^^i^^'^^i^^'^'^ utilized to construct 
high-fidelity quantum gates from a set of classical con- 
trols, which may be pulsed in counter-intuitive ways, that 
would be difficult to find with traditional analytical tech- 
niques. 

There are many different ways in which to formulate 
the optimal control problem, in terms of the target quan- 
tum operation, and in terms of numerical search meth- 
ods. In a very broad way, the general idea of optimal 
control methods is to iteratively search for a target set of 
control field parameters that drive a system Hamiltonian 
in a controlled manner. Depending on the system con- 
sidered, the target parameter could be quantum process 
fidelity, or operational time of the evolution, or the value 
of a specific quantum observable. 

One such example of an optimal control method is the 
standard numerical algorithm known as gradient ascent 
pulse engineering (GRAPE) optimal quantum gates 
(in terms of process fidelity) are constructed by itera- 
tively varying control parameters, and analytically calcu- 
lating the gradient changes, in order to find the minimum 
fidelity points in the underlying unitary space. These 
techniques have been applied to ion traps^^ vibrational 
modes in molecules;^ as well as in the superconducting 
].ggjj^Q56j58 Tff[f^]j significant success. 

Our future work will involve utilizing these numerical 
techniques to find smooth control field parameters that 
modify the above simulated BSB pulse to one that elim- 
inates leakage to higher order field modes. Hopefully, 
once we eliminate the final issue of qubit leakage, we 
will utilize these numerical algorithms in full open sys- 
tem calculations to develop high fidelity control pulses 



for superconducting systems undergoing decoherence. 



VII. CONCLUSIONS 

Our simulations have confirmed a usable control, tar- 
get, and BSB pulse matrix, as well as the required 
operational pulse rise/fall time envelope published pre- 
viouslyf^ for an experimentally successful microwave 
pulse-induced C-NOT gate. However they have also 
shown that this set of pulses, in particular the BSB set, 
incites unwanted qubit leakage, at a minimum, into the 
first and second excited states of the coupling resonator, 
which makes them troublesome for usage in a more QIP- 
focuscd architecture. 

Our proposed solution to this problem, replacing the 
BSB pulse sequence with a sequence based on optimal 
quantum control concepts, seems to have the potential to 
reduce the observed qubit leakage. Yet it is still uncer- 
tain whether or not the required pulse sequences needed 
for these codes can be realized experimentally. The main 
reason is that our current architecture has a 1:1 reso- 
nance between qubit number and microwave lines. This 
was done in order to avoid cross-talk between qubits. 
More complicated pulse sequences on this architecture, 
such as those needed for GRAPE codes, run the risk of 
initiating cross-talk between qubits and other microwave 
lines, or exciting on-chip fluctuators close to, or within, 
each qubit's junction. 

Currently, we are planning to explore these issue in 
more detail by investigate other GRAPE-inspired multi- 
qubit /resonator gate sequences, as well as modifications 
to the architecture. Future work will also look into more 
accurately modeling the dynamical impact of the various 
low- and high- frequency noise sources inherent in our 
system;^ using a more generalized open quantum system 
treatment 
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